Radiative corrections to the background 
of jjL — > decay 



A.B. Arbuzov 1 , O. Krehl 2 , E.A. Kuraev 1 , 
E. Magar 1 , and B.G. Shaikhatdenov 1 
1 Joint Institute for Nuclear Research, Dubna, 141980 Russia 
2 Institut fur Kernphysik, Forschungzentrum Julich GmbH 
52425 Julich, Germany 

Abstract 

Radiative muon decay in the kinematics similar to the neutrinoless decay 
u — > e7 is considered. Radiative corrections due to one-loop virtual pho- 
tons and emission of additional soft or hard photons are taken into account. 
Analytical expressions and numerical estimations are presented. 
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1 Introduction 



Since the discovery of the muon in 1936 its relation to electron is a puzzle. Really, 
the only difference between these two elementary particles is in their masses. The 
lepton number conservation law has no deep sources in space-time properties or 
gauge theories. Moreover many extensions of the Standard Model predict processes 
with violation of this law (fi — > ej, ej'j, eee etc.). Intensive search of these 
extensions was performed in 1977 0. The modern state of the subject is elucidated 
in papers |3| and references therein. Indeed, if there is a unification of quarks 
and leptons, then the existence of b — > 37 decay leads to that of p, — > cy. Different 
models give a wide range of predictions for the branching ratio of this neutrino less 
muon decay. The present experimental upper limit || on the branching ratio is 



fi = ^-^<4.9.10- 11 . (1) 
riot 



This value imposed already strong restrictions on parameters of supersymmet- 
ric pi |5| and other models |]]. In the model independent approach J7| one gets 
boundaries on parameters of possible structures in the matrix element of the muon 
decay. Several new experiments are planned to improve the precision. They will 
either find the decay or put much more stronger restrictions and even discriminate 
some models. The forthcoming experiment at PSI (if doesn't find the decay) will 
put the limit on the // — > e 7 decay branching ratio of about 5 • 10~ 14 . Another 
experiment is proposed at BNL, where they are going to reach the level of 10~ 16 . 
These experiments are very important, since they have rather wide possibilities for 
the search of new physics comparable with those of high energy colliders. In this 
paper we consider the important background process 

M -> e{p 2 ) + 7 (fcx) + (z/ M + v e )(q) (2) 

in the kinematical situation, imitating the neutrinoless decay. Namely, we suppose 
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where q is the 4-momentum carried by neutrinos, and M is the muon mass. The 
width in the lowest order of perturbation theory was calculated many years ago ||. 
The expression for the width reads: 
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q = p -P2 - h, ^i = k°. (3) 



The validity of this formula may be confirmed in the limiting case of soft photon. 
A multiplier 2 was lost on right hand side (rhs) of expression for the width in ||. 
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The polarized muon radiative decay was considered in || , and as a background to 
the neutrinoless decay it was extensively discussed in Ref. 0. 



2 Radiative corrections 



In the imitating kinematics (IK) we introduce the relative energy deviations of the 
hard electron and photon from M/2 and the acollinearity angle 9: 
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and expanding the expression for the width in the Born approximation 
obtain: 
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Here s denotes the spin of the muon, and (p is the azimuthal angle between planes 
formed by (s,p 2 ) and (s,ki) in the rest reference frame of the muon. Note that 
averaging the above expression over the angle ip leads immediately to the result 
presented in ||. We shall name higher than second order contributions on the rhs 
of (|3p (and Si in rhs of (|J)) as relativistic corrections. In this paper we will consider 
the radiative corrections to this width bearing in mind virtual corrections described 
by the Feynman diagrams drawn in Fig. 1 together with those arising from the 
emission of additional soft and hard photons. 

For the measurement of an additional hard photon emission, two cases have to be 
considered: with and without external magnetic filed. In the case without magnetic 
field the additional hard photon, moving along the final electron trajectory within a 
small angle, which is equal to the detector angular resolution, is registered together 
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Figure 1: The subset of Feynman diagrams for radiative muon decay. 



with the electron. In the opposite case (with magnetic field) those events will be 
rejected from statistics due to criterion: the energy of the electron is less than 
the maximum energy within some accuracy. The standard calculation of one-loop 
virtual corrections can be considerably simplified by using the IK features. Some 
details of our calculations (traces, vertices and the Tables of relevant integrals) are 
given in Appendices. 

Ultraviolet divergences of loop integrals are eliminated in a standard way using 
the renormalization constants of the wave functions of electron and muon: 
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where m, A, A are the electron mass, infrared and ultraviolet cut-off momentum 
parameters, (K<m,A» M). The final result for the one-loop virtual corrections 
reads: 
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Taking into account the emission of additional soft photon requires some care. The 
reason is that the energy-momentum carried by soft photons as well as by neutrinos 
cannot in principle be distinguished in the experiment. We introduce some small 
energy fraction parameter Ai = 2ui soit /M ai,a 2 ,9 which should not affect on 
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observable quantities and actually cancels out in the final result. Emission of an 
additional soft photon, having energy lesser than MAx/2, can be taken into account 
in a usual way [|K|. The corresponding expression looks as follows: 
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Let us suppose, that a photon with momentum k 2 , having energy more than w so ft, is 
emitted in such a way that we have still allowed values of the final electron and hard 
photon momenta. In this case the additional photon cannot be called soft, because 
it changes the kinematics of the process. We have to consider the corresponding 
contribution applying complete set of kinematical restrictions. The main condition 
is that the missing momentum squared must be positive: 



(p -P2 ~ h) 2 > 0, q = q + k 2 . 



(10) 



Having in mind that the matrix element squared is proportional to the second power 
of small neutrino momenta, we can write down the contribution under consideration 
in the factorized form 
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The difference in respect to the case of pure soft photon emission is that we have 
the shifted quantity R instead of the Born one (R) under the integral sign. The 
above expression guarantees that the energies and angles of the observed electron 
and photons are the same as defined in (H) . Transforming the above formula we get 
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In this expression we introduced an auxiliary parameter 9q in order to separate 
the contribution, when the additional photon is emitted collinear to the electron 
momentum; S C 1. So, the first term of Eq. fljjD can be integrated analytically 
in order to keep track of the leading logarithmic part. We checked that the final 
expression does not depend on 9 . 

Then we arrive to the total answer, that has the form 

dr 

— = 1 + 5 1 + 5 V + 5 S + 5 7 • (13) 

Cll o 

The dependence on the soft photon parameter Ai cancels out in the sum 5s + 5 y 
whereas the fictitious photon mass A disappears in the sum 5s + <V- 

If the experimental set-up does not distinguish in the detector an electron with 
a collinear photon, we have to modify our results in the following way. Let 9 define 
the aperture of the narrow cone, within which the two particles would be detected as 
a unique one. Then we should take the non-shifted value for R in the first integral 
of Eq. ([12]). We have to add also the rest contribution of hard photon emission 
within the same cone. It can be obtained using the quasireal electron method [|TT|j : 

dr- dr * / dl i + (i-*)\(MM (14) 



daido^d^ da 1 da 2 9d9 n J x \ 2m 

•''max 

*W " — • 

The lower limit comes here from the B-function in the first integral of Eq. flT2[) . 

In the presence of a magnetic field, when the electron trajectory is curve, the 
above expression will give a part of the background to the process fi — > e^yy, con- 
sidered in paper [I2|. Really, the final electron will have the energy M(l — x)/2, 



whereas the quantities <j\ , 02 , 9 are the same as in the case of single photon emis- 
sion. 



3 Conclusions 

In Table 1 we give numerical values for 5±, 5sv-y = 5s + 5y + 5^ versus 01, 02, 9, 
which characterize the missing energy and momentum (see Eq. (f|)), and the degree 
of muon polarization £. For typical expected values of a\ ~ a 2 ~ 9 ~ 1CT 2 (we give 
6 points) one can see, that the relativistic and QED corrections should be taken 
into account on the same footing. The resulting correction to the Born-level decay 
width dr (see Eq. (]6|)) is given by the sum 5sv^ + <5i- 

A measurement of the radiative muon decay in the kinematics close to that of 
neutrinoless decay is required to get an independent normalization for the search 
of the decay /i — > ej. For this aim our results, we hope, would be important. 
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4.9 


-13.1 


-13.2 
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4.0 
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-6.5 
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-7.4 



Table 1: Numerical estimations for the corrections 5\ and dsv-y versus a±,a2,9,^ 



We would like to mention here result obtained in on the background to the 
three lepton neutrinoless decay [i + — > e + e + e~. For an experimental set-up when 
the electron and positron energies are measured, it reads 

.( 2 - w ) 2 ln— , «,= (£+ + £+ + £), «w 2. (15) 



T dA tt 2 36 v ' m 2 ' M 

We have to discuss some features of the results presented. At first we note, that 
the large logarithm L does not factorize before the Born-like structure (R), as one 
may expect. We claim that the factorization theorem, which was proved for high 
energy processes, should not work here. Another problem is that if one integrated 
out over the whole phase volume of the second photon, he would still have in the 
answer the logarithm of the mass ratio. Formally, this violates the Kinoshita-Lee- 
Nauenberg theorem fI3|; the formula is infinite in the limit m — > 0. But again, the 



conditions of the theorem allow us to say, that the process of radiative muon decay 
is a legal exception. One can see the same situation in radiative muon decay at the 



Born level [fLOL 15 
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Appendix A. Tables of integrals 

Here we put the tables of relevant integrals appearing in the loop momentum inte- 
gration. The denominators of amplitudes, which correspond to Feynman diagrams 
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2 ..,2 ... 1.2 07.., , 1 72 (A.l) 



drawn in Fig.l, have the following form: 

(1) = (P-k) 2 -M 2 , (2) = (p 2 - k) 2 - m 2 

(1) = (p — fci — A;) 2 - M 2 « A; 2 - 2A;p 2 - M 2 , 

(2) = + fci - A;) 2 - m 2 « A; 2 - 2A;p + M 2 , 

(0) = A: 2 - A 2 . 

We use a symbol ~ to underline the peculiarity of imitating kinematics. Namely, 
working out traces we use 

p\ = k\ = 0, 2p 2 A;i = M 2 = 1, q = 0. (A.2) 

The scalar integrals considered have a form 

dk r dk r dk „ d 4 A; 



dA; = — . (A.3) 



mr J mmr j (mmr 

Vector and tensor integrals are parametrized as follows: 

/ = ck ^ + d?2 ' 

r k^k u dk 

j = g<r + aK k + + i(h, P 2r , 

J k^Wdk = ^ + G (2) ( ^ ft) + K(A;i) 3 + r(p2 )3 + 

V(P2,^l^l)+p(P2,P2^l)r CT , (A.4) 

where we denote different symmetrical combinations, for instance: 

(0, a p = ^a* + + , (a, = a l V + a J 6\ ... (A.5) 

Below we put the values of the coefficients and the scalar integrals. In the tables 
2 + 7 we used Y = In ffc, L = In M ? x = f , Z = In M™. All the integrals we put 
in dimensionless form by setting M — 1. 



Appendix B. Gauge invariant subset of Feynman 
diagrams 

Amplitudes, describing Feynman diagrams with loop correction to the real photon 
emission vertex, and the ones, taking into account self- energy of fermions (typical 
diagrams are shown in Fig.la,b), provide a gauge invariance in respect to the real 
photon polarization vector. It has a universal form and may be taken into account 
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by substitutions in Born amplitude of the form 
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Li 2 (l)-Li 2 (^) , Li 2 (,) = -/ o 



nl ln(l — zx) 



x 



dx, 



where e is the polarization vector of the real photon. In the IK we have A\ = 
1/(4M), A 2 = (n 2 /Q — l/4)/M 2 . In a similar fashion for the diagrams, which can 
be obtained from depicted in Fig.la,b by emitting a real photon from another leg, 
we have 



_ . , p 2 + k\ + m a _ . . 
U{P2) 2p 2kl 



B 1 (e-fci-^-J +B 2 he 

V P2K1/ 



(B.2) 



In the IK, omitting the terms disappearing in the limit of zero electron mass, we 
have B 1 = 0, B 2 = (1/2 - 2L)/M 2 . 



Appendix C. Averaging on neutrino states, traces 

To rearrange bispinors in the matrix element we use Fierz identity: 

uiO a u 2 u 3 O a u 4 = -u 3 O b u 2 UiO b U4, O a = 7 (1 + 7 5 )/2. (C.l) 
Summing on the neutrino spin states of the matrix element squared one obtains 

Eu 3 O a u 2 (u 3 O b u 2 )* = 2(( qi q 2 ) ab - qi q 2 g ab ) = L ab . (C.2) 



/ 



Averaging over the neutrino momentum is performed using the invariant integra- 
tion: 

d ^^ 1 #(gi + q 2 - q) = l{q 2 g ab + 2gY). 
<?10<?20 6 

Application of this formula to the tensor L ab gives the result: 

l3„ i3 



(C.3) 
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The doubled interference term of Born and one-loop amplitudes looks as follows 
(we consider IK): 



•_>>:. u;,.u, 



2 8 aG 2 F n 
3M 



—A{T\ - A 2 T 2 + B 2 T 3 + 



5/* 



Si S 2 S3 S4 



.(0)(2)(1) (0)(1)(2) (0)(1)(2)(1) (0)(1)(2)(2) 
where the traces are: 



(C.5) 



Tj = ^Tr{f t ac )O ac , Si = -Tr(S% - S%)O ac . (C.6) 

The list of T™, Sg, S% is given below: 

fi = p 2la { lb -2k m /M 2 ){p + M) lb {p 2 + M) lc /M\ 
f 2 = P + 2 la kMP + M) lb {p 2 + M) lc /M\ 

% = p2lakllbPlcPlcPlb/M i , 

S31 = (p + l)7cP7bP>)v(P2 - k)laip2 - k + l)7b(p - k + 1)7^ , 

^12 = (P2 + l)7fc(P+ l)7cP7feP27^(P2 - A;)7a(P2 - fc + 1)7^ , 

5 21 = (p + lhcPlbp2lbPlv(P - k)la(P ~k+ 1)7^ , 

S41 = (p + 1)7cP76P27^(P2 - A;)7b(p - fc)7a(p -k + 1)7^ > 

£32 = (p+ 1)76 (P2 + l)7cP27 M (P2 - k)7a(P2 - fc + l)7fc(P - k + 1)7m> 

#42 = (p+lhb(p2 + ^hcp2l^(p2-k)^ b (p-k)-f a (p-k + l)-f fl , 

Su = (p 2 + l)7b(p + l)7b(P2 + 1)7^27^2 - k)-i a (p 2 -k + 1)7^ , 

522 = (p+ 1)76 (P2 + ^)lcp2lbPl l ,{P - k)laip - k + 1)7^, 



5*1 — 5*11 — 5*12, S 2 — S 2 i — S 22 , S3 — S31 — 532, S*4 — 5*41 — S42. (C.7) 
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Table 2: Scalar integrals with 2 denominators 
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Table 3: Scalar integrals with 3 and 4 denominators 
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Table 4: Vector integrals with 2 denominators 
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Table 5: Vector integrals with 3 and 4 denominators 
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Table 6: 2-rank tensor integrals with 3 and 4 denominators 
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Table 7: 3-rank tensor integrals with 4 denominators 
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